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Abstract 

Strongly confining an ultracold atomic gas in one direction to create a quasi- 
2D system alters the scattering properties of this gas. We investigate the 
effects of confinement on Feshbach scattering resonances and show that strong 
confinement results in a shift in the position of the Feshbach resonance as a 
function of the magnetic field. This shift, as well as the change of the width of 
the resonance, are computed. We find that the resonance is strongly damped 
in the thermal gas, but in the condensate the resonance remains sharp due 
to many-body effects. We introduce a 2D model system, suited for the study 
of resonant superfluidity, and having the same scattering properties as the 
tightly confined real system near a Feshbach resonance. Exact relations are 
derived between measurable quantities and the model parameters. 
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I. INTRODUCTION 



The ability to tune the interatomic interaction strength by using a Feshbach resonance 
has led to renewed interest in the scattering theory of cold trapped atoms, in particular 
in relation with the realization of resonant superfluidity of fermionic atoms [1]. Also the 
ability to confine Bose-Einstein condensates in periodic potentials has attracted much in- 
terest recently, since the resulting optical lattices are ideally suited to probe fundamental 
quantum effects such as the superfluid to Mott insulator transition [2]. One-dimensional 
optical lattices allow to create stacks of 'pancake'-shaped clouds filled with neutral atoms 
[3-5], where an effectively two-dimensional (2D) regime can be achieved. In this Letter, we 
address a central question on the crossroads between the two aforementioned achievements, 
Feshbach scattering and optical lattices: what happens to the Feshbach resonance as the 
atoms are strongly confined and become (quasi-) two-dimensional ? 

In Sec. II, we derive the 2D resonant scattering amplitude in terms of quantities that 
can be measured in experiments on three-dimensional (3D) gases. For many-body calcula- 
tions, it can be helpful to replace the complex real atomic interaction by a simpler effective 
interaction. To achieve this, we investigate in Sec. Ill a schematic model of coupled atoms 
and molecules that interact through a contact potential. Such model systems already have 
been used before in 3D [6], but the connection between the model parameters and the ex- 
perimental quantities in 2D is not yet available. By matching the scattering properties of 
this model to the results from Sec. II, we derive expressions for the strengths of the effective 
contact interactions as a function of experimentally measurable (and tunable) quantities. 
This solves the problem of determining the parameters of a Friedberg-Lee like Hamiltonian 
[7] for atomic gases. Such models have been studied in the context of high temperature 
superconductivity, but for these systems fixing the model parameters remained a problem. 

II. QUASI 2D RESONANT SCATTERING 

The existence of a Feshbach resonance relies on the presence of a bound state in a 
closed channel, with an energy equal to the energy of the colliding particles. In the case 
of trapped cold atoms, the closed channels involve states with higher Zeeman energy in the 
external magnetic field. The energy of the bound state (molecule) is tunable by changing 
this magnetic field. Because these molecular states are small compared to the length scale of 
the confinement potential, they are not expected to be influenced strongly by the trapping 
potential. However, this is only a simplified picture. Because of the coupling between the 
molecular state and the open channels, the binding energy of the molecule acquires a complex 
shift (a change of the binding energy and a decay width). This shift depends on the nature 
of the open channels and consequently on the tight confinement. 

If the external confinement potential varies slowly on the scale of the atomic interaction 
range, the two-body scattering can be treated analytically, even in the vicinity of a Feshbach 
resonance (see also Ref. [8]). In the following, we denote the range of the van der Waals 
interatomic potential with R e and the oscillator length of the tight confinement in the axial 

direction by £ z = ^Jh/ (2m Te( iUJ z ), with m re d = {jn^ 1 + m 2 *) the reduced mass of the 
scattering atoms and u z the characteristic frequency of the trap in the axial direction. For 
interatomic distances R e <C r <C £ z , both the interatomic and the confinement potential are 
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negligible compared to the total energy and the asymptotic s-wave scattering wave function 
is consequently that for free 3D motion 
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sin (qr + 5 3D (E)) 
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(1) 



with q = y/2m TC dE /h. The energy E is measured in the center of mass frame of the scattering 
atoms, relative to the energy of two atoms in the open channel at infinity. The phase shift 
Ssd (E) is determined by the solution of the free 3D scattering problem. At a Feshbach 
resonance, this phase shift shows resonant behavior, but the wave function (1) remains 
valid. Close to the resonance, the energy dependence of S 3 £> (E) is important, because it 
characterizes the resonance. From the analytic theory of Feshbach resonances [9] it follows 
that it is given by 



Ssd (E) = 5b g — arctan 
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where 5 bg is the 3D 'background' scattering phase shift related only to the open channels, and 
Tf,Ef and A^ are respectively the width and the position of the 3D Feshbach resonance 
and the shift of that position due to the coupling with the open scattering channels. 

Next, we look at large interatomic distances r » E e , where the motion is governed by 
the confining potential and must be a solution of 
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Petrov and Shlyapnikov have analyzed this problem [10]. Using their solution of Equation 
(3) and matching it for r — > to (1) we find for the scattering amplitude 
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where 5 3D (E) is given by (2) and C = 0.915 [10]. Making use of formulas derived in [9], 
cot Ssd {E) can be related to measurable quantities 



cotS 3D (E) = 



1 E-{B-B Q + AB{qa bg f}A^ 
' qa bg E — (B — B — AB) A/x : 



(5) 



where B and AB are respectively the 3D experimental position and width of the Feshbach 
resonance, and a^ g is the background scattering length. A/x = /Xi + /x 2 — /x res ; /ii i2 are the 
atomic magnetic moments and /x res is the magnetic moment of the resonant state. 

Formula (5), together with (4) gives the 2D resonant scattering amplitude for a quasi 2D 
situation in terms of parameters that can be measured in 3D experiments. It should be noted 
that for low energies (small compared to Ep + A^, which is typically of the order /x# x 1 G 
= ksx 67.17 /iK, where /x# is the Bohr magneton), the energy dependence of (5) is weak and 
cot 5 S d (E) is well approximated by its value for E — 0, namely — q£ z cot S SD (E) — > £ z /a 3D , 
which would yield the usual relation a^ D = a bg [l — AB/(B — B )] for a Feshbach resonance 
in 3D. 
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A remarkable feature of the expression (4) is that it shows that in two dimensions the 
scattering amplitude never diverges, because \f 2 u {B)\ 2 < 16, irrespective of E or Ep. This is 
in contrast with the 3D case where the imaginary part in the denominator of the scattering 
amplitude (= T F ) is proportional with y/~E and vanishes as E — > [9]. However, when 
many-body effects are taken into account for bosonic atoms, the damped character of the 
Feshbach resonance in two dimensions disappears. Lee et al. show in Ref. [11] that in the 
limit of zero temperature in a BEC the relevant matrix element of the 2D many-body T- 
matrix Tmb is to a good approximation equal to the analytic continuation of the two-body 
T-matrix T 2b (E) = (h 2 /m) f 2D (E), given by (0| T MB (E = 0) |0> = (0| T 2b (-A) |0> , where 
A is the chemical potential of the Bose-gas. By this procedure, the in term in expression 
(4) disappears, so that a sharp resonance is expected in a Bose-Einstein condensed gas at 
low temperatures. Formulas (4) and (5), imply that the resonance position in the quasi-2D 
condensate, Bq D , is different from B . The shift Bq D — B depends on the chemical potential 
and diverges for A = e* (uj z ) = Chcd z /n exp[\^2irl z /a b g], which is the energy of the weakly 
bound 2D state if only the background scattering is taken into account [10]. At this value 
of the scattering energy, also the width of the resonance in a BEC diverges. For the thermal 
gas, there are no coherent many-body effects and \f 2 p> {B)\ 2 (expression (4)) is the relevant 
quantity, for example in the calculation of the rethermalization rates. 

The scattering amplitude in the neighborhood of the Feshbach resonance at -Bo = 154.9 
G in a 85 Rb condensate (experimentally explored in Ref. [12]) is shown in Fig. 1 as a 
function of the confinement frequency and the magnetic field. The resonance parameters 
are AB = 11.9 G, a bg = — 380ao? A/i = — 2.23/i# [12,13]. The chemical potential was fixed 
at \/h = 5 kHz. The white regions in Fig. 1 correspond to |/(— A)| > 100 and indicate 
the location of the resonance. The divergence of the shift at e* (uj z ) = A can clearly be 
seen. Inset (b) of Fig. 1 shows the scattering amplitude f 2 o{— A) for condensate atoms 
at uj z = 0.4 MHz, as a function of magnetic field. For this value of u z , we find a shift of 
Bq D — Bq = —16.5 G, indicated by an arrow in inset (b). 

In experiments, it is the radial confinement and the number of particles that are con- 
trolled, rather than the chemical potential. The chemical potential for a homogeneous system 
is equal to A = (0| T M b (E = 0) |0) x n b , with n b the density of bosons and accordingly has to 
be determined self-consistently. When the chemical potential is determined self-consistently 
in the Thomas- Fermi approximation, again we find that the resonance is shifted. But to the 
left of the resonance no self-consistent solution for A can be found. In this region we expect 
the interactions to become attractive, and the absence of a self-consistent solution can be 
linked to the instability occurring in condensates with attractive interactions. 

In Fig. 2, |/ 2 d {E)\ 2 , the squared modulus of the scattering amplitude for the thermal 
cloud, is plotted for E/ks = 100 nK. For the thermal gas the maximum value of the 
scattering amplitude is bounded. 

III. 2D CONTACT SCATTERING MODEL 

To avoid having to take into account the details of the real interaction, it is convenient 
to introduce a 2D model with an effective contact potential that has the same scattering 
properties as the real interatomic potential in the quasi-2D system. This effective potential 
can then be used to tackle many-body problems such as the superfluidity of a two-component 
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fermion gas. Kokkelmans et al. have performed such a study for the 3D case [6]. They 
introduced a molecular degree of freedom and showed that the model of atoms and molecules 
that are converted into each other leads to a scattering length of the Feshbach type. We 
construct an analogous model to describe resonant scattering in a quasi-2D gas. Our model 
2D Hamiltonian for the scattering of two different particles 1 and 2 (bosons or fermions) 
takes the form 
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(6) 



with z/ the binding energy of the molecule and g m is the parameter related to the formation 
of a molecule out of two atoms. V p is the interaction far from resonance and ipi,ip2 and 
■ipm annihilate respectively the atoms 1, 2 and molecules. We solve the scattering problem 
for two atoms in the center of mass frame along the lines of Ref. [6] and we find for the 
scattering amplitude 



f 2D (k,p) =2V2tt 



2 v / 27( 



n 2 k 2 

2m. r , 



V 



gmred y P ( ?i 2 k 2 
h 2 \2m Icd 



9n 



1 2 2m. 



+ 



2tt 



In 



k 2 



- 1 



+ m 



(7) 



which is independent of p. K c is a momentum cutoff, introduced to keep the integrals over 
p finite. Expression (7) has the same structure as (4). Comparing both formulas, using 
(5), gives the relation between the physical quantities and the model parameters. With the 
notation x — 2A_BA/ia 2 „m rc d/h 2 , we find 



V F = 



1 \P2j\h 1 dbg 



1 - x m Icd 4 

for the non-resonant interaction. The molecular binding energy is 

v=(B-B )A»/(l-x) 
and the coupling strength is given by 

\g m \ 2 = V p (ABA^ + xv)- 



(8) 



(9) 



(10) 



From the comparison of (4) and (7), it follows that the cutoff momentum K c is related to 
the tight confinement harmonic oscillator length by K c - 



/C/(ir£ 2 z ). The right hand side 
of expression (10) can be negative for some values of x an d B. Because the description in 
terms of atoms and molecules describes the physics close to the resonance, we need certainly 
that for B = Bq the RHS of (10) is positive. From the theory of reference [9] it follows that 
a^gABAfi > 0, so our model (6) is applicable only for % < 1. For the case of 85 Rb around 
B = 154.9 G, we find that x — —126.3 < 1 so that the model is applicable. In contrast 
with the 3D case, the parameters V p and g m in 2D do not depend on the cutoff momentum 
and need no renormalization, but the cutoff momentum itself is fixed and depends on the 
tight confinement frequency. Thus, formulas (8), (9) and (10), together with the knowledge 
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of the experimental quantities a^, Bq, AB and A/x from the 3D Feshbach resonance, fix 
unambiguously the parameters of the model Hamiltonian (6). This allows us to describe 
a quasi-2D mixture of gases that are tightly confined in the axial direction by a harmonic 
potential with characteristic length £ z . As a first check of the model we have found that far 
from resonance this contact model with a finite cut off reproduces the critical temperature 
for superfluidity derived in Ref . [14] . 

IV. CONCLUSIONS 

The scattering of cold neutral atoms that are strongly confined in one spatial direction is 
studied analytically close to a Feshbach resonance. We find that there exists a confinement- 
induced shift in the position of the Feshbach resonance as a function of the magnetic field, 
and show that for a 85 Rb condensate this shift is experimentally detectable. We also find 
that the Feshbach resonance in the confined thermal gas is damped, in the sense that the 
scattering amplitude cannot diverge. Based on the present analysis of Feshbach scattering 
in a quasi-2D gas, we set up a Friedberg-Lee like model system with an effective contact 
interaction, suitable for the study of resonant superfluidity in confined gases. Exact relations 
are derived linking the model parameters to experimentally measurable quantities. 
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W.O.G. project WO.025.99N.and the GOA BOF UA 2000, IUAP. 
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FIGURES 



FIG. 1. The scattering amplitude fiD (— A) for a quasi-2D 85 Rb condensate, with chemical 
potential A/ft = 5 kHz, is shown color-coded as a function of the magnetic field and the trapping 
frequency in the tightly confined direction u z , close to the 3D Feshbach resonance at Bq = 154.9 
G. The white regions correspond to |/(— A)| > 100 and contain the location of the resonance in 
the quasi-2D condensate. In the insets (a) and (b) Jid (— A) is plotted at fixed magnetic field and 
fixed frequency, respectively, as indicated by the lines (a) and (b). At a confinement frequency of 
uo z = 0.4 MHz (inset (b)), the Feshbach resonance is shifted by —16.5 G, as shown by the arrow. 

FIG. 2. The modulus square of the scattering amplitude \f(E)\ 2 for a quasi-2D 85 Rb thermal 
cloud is shown for T = 100 nK close to the Feshbach resonance at Bq = 154.9 G for and to z = 10 5 
Hz. Results are shown with (full line) and without (dashed line) taking thermal broadening into 
account. In both cases, the divergence of the scattering amplitude has disappeared, due to the 2D 
character of the gas. 
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FIGURE 2 



